y=b/sqrt(1-(a/sqrt(1-y^2/2))^2/2)